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Motivated by recent experiments on an S = 1/2 antiferromagnet on the kagome lattice, we inves- 
tigate the Heisenberg Ji — J2 model with ferromagnetic Ji and antiferromagnetic J2. 
Classically the ground state displays Neel long- range order with 12 noncoplanar sublattices. The 
order parameter has the symmetry of a cuboctahedron, it fully breaks SO(3) as well as the spin flip 
symmetry, and we expect from the latter a Z2 symmetry breaking pattern. As might be expected 
from the Mermin- Wagner theorem in two dimensions, the SO(3) symmetry is restored by thermal 
fluctuations while the Z2 symmetry breaking persists up to a finite temperature. 
A complete study of S — 1/2 exact spectra reveals that the classical order subsists for quantum 
spins in a finite range of parameters. First- order spin wave calculations give the range of existence 
of this phase and the renormalizations at T — of the order parameters associated to both symme- 
try breakings. This phase is destroyed by quantum fluctuations for a small but finite J2/IJ1 — 3, 
consistently with exact spectra studies, which indicate a gapped phase. 



THEORETICAL AND EXPERIMENTAL 
ISSUES 



Whatever the nature of the spin, classical or quantum, 
the first neighbor Heisenberg antiferromagnet on the 
kagome lattice fails to display Neel-like long-range 
order. Classically, it is characterized by an extensive 
entropy 1,2 at T = 0. Quantum mechanically the spin-1/2 
system has an exceptionally large density of low lying 
excitations 3 ' 4 reminiscent of the classical extensive 
entropy. It is still debated whether and eventually how 
this degeneracy is lifted in the quantum limit 5,6 . 
An essential issue concerns the influence of pertur- 
bations: classically the effect of a second neighbor 
coupling J2 has been very early studied by Harris and 
co-workers 7 . They showed that an infinitesimal Ji is 
sufficient to drive the system toward an ordered phase 
with the three spins around a triangle pointing 120° 
from each other. Antiferromagnetic second-neighbor 
coupling ( J2 > 0) favors the q = Neel order of this 
pattern on the Bravais lattice, whereas there are nine 
spins per unit cell for J 2 < (q = \fi x v3 order). The 
effect of Dzyaloshinsky-Moriya interactions has also 
been analyzed 8 . To our knowledge the reduction of the 
order parameter by quantum fluctuations has only been 
studied through exact diagonalizations 9 . This approach 
points to an immediate transition from the "disordered 
phase" at the pure J\ > point, to the semiclassical 
Neel phases. 

Up until now the J\ — Ji model on the kagome lattice 
has only been studied for antiferromagnetic J±. Many 
magnetic compounds 10-13 with this geometry have been 
studied so far, but most of them have spin 5 = 3/2. A 
few compounds with S — 1/2 Cu ions have recently been 



synthetized 14-16 . None of them can be described by a 
pure isotropic first neighbor antiferromagnetic Heisen- 
berg model. Recent experimental work on an organic 
compound with copper ions on a kagome lattice 1 gives 
indication of competing ferromagnetic and antiferromag- 
netic interactions. 

It is thus the purpose of the present work to extend 
the previous study of the J± — Ji model to ferromag- 
netic nearest neighbor coupling (Ji < 0). The Hamilto- 
nian reads as 



H — Ji y ' Si • Sj + J2 ^ ' Si ■ S 



A'- 



(1) 



<i,j> 



«i,k» 



where the first and second sums run, respectively, over 
pairs of nearest neighbors < i,j > and next-nearest 
neighbors << i, k >>. 

For a pure ferromagnetic J\ coupling the system is indeed 
in a ferromagnetic phase. For a pure antiferromagnetic J 2 
interaction the model reduces to three decoupled kagome 
lattices with antiferromagnetic interactions and has thus 
an extensive entropy in the classical limit. The behavior 
of the model between these two limits is the object of the 
present study. 

In Sec. II the classical ground state of the Hamiltonian (1) 
is investigated and the phase diagram of the model is 
given in the classical limit. For competing interactions 
Ji < and Ji > the model exhibits an ordered phase 
with 12 sublattices that fully breaks SO(3) as well as a 
discrete symmetry (chiral symmetry breaking). We show 
that contrary to the Neel order, which breaks a contin- 
uous symmetry and therefore is destroyed down to in- 
finitesimal temperatures, the chiral order survives ther- 
mal fluctuations and undergoes a phase transition at fi- 
nite temperature. 

In Sec. Ill we study the S — 1/2 quantum model using 



exact diagonalizations and show the premise of the semi- 
classical ordering on samples up to 36 spins. 
In Sec. IV we study the effect of long wavelength quan- 
tum fluctuations on this semi-classical order in the spin- 
wave approximation. It appears that the twelve sublattice 
Neel order survives quantum fluctuations in an extended 
range of parameters but is destroyed for a non zero J2 > 
(■Wi|>3). 

In Sec. V we show that exact diagonalizations in this 
range of parameters indeed point to (a) gapped phase (s). 



II. CLASSICAL APPROACH 

In this section we restrict ourselves to classical spins: 
the spins Si are usual unit vectors living in a three- 
dimensional space. 



A. Ground state for Ji < and J 2 > | Ji | / 3 

To investigate the nature of the ground state of the 
Hamiltonian (1), we first Fourier transform it to find its 
lowest-energy modes. The kagome lattice having three 
sites per Bravais cell, we get three branches. For J\ < 
and J2 < we find a single minimum at q — corre- 
sponding to the expected ferromagnetic ground state. 
Upon increasing J 2 > 0, we find three degenerate min- 
ima at the edge centers of the first Brillouin zone, the 
classical transition occuring at J% = |Ji|/3. The three 
modes q = Xi, 2 ,3 (Fig. 1) are the only solution as long 
as J 2 > \Ji |/3 and in the limit of pure J2 > one recov- 
ers the flat zero-energy branch of the pure J\ > case: 
we then have three decoupled kagome lattices with near- 
est neighbor coupling J 2 . In the J\ < 0, J2 > |Ji|/3 re- 





FIG. 2: (Color online) Classical order parameter for Ji < 
and J2 > |Ji|/3 with the associated lattice symmetry break- 
ing. The 12 sublattices and the 4 different types of hexagons 
are numbered. Also shown is the plane containing the six spins 
of hexagon 1 (dashed) . 



their nearest neighbors. We thus have four hexagons in 
the unit cell defining four different planes oriented like 
the faces of a tetrahedron. The polyedron whose vertices 
coincide with the directions of the sublattices is named a 
cuboctahedron, and we will refer to that order parameter 
as the cuboc phase in the following. Consistently, Monte 
Carlo simulations with a Metropolis algorithm reveal a 
local cuboc Neel order with fluctuations increasing with 
temperature. 

We now have a complete picture of the classical phase 
diagram of the Hamiltonian (1) at T — in the entire 
Ji — Ji plane (Fig. 3). 




FIG. 1: (Color online) The first Brillouin zone of the kagome 
lattice with its points of high symmetry: the center of zone 
r (q = 0), the two zone corners Mi and M2, and the three 
edge centers Xi, X2, and X3. 



FIG. 3: (Color online) Classical phase diagram of the Heisen- 
berg Ji — J2 model on the kagome lattice at T = 0. The tran- 
sition between the ferromagnetic state and the antiferromag- 
netic cuboc state occurs at J2 = — Ji/3. 



gion, the unit cell compatible with the three edge centers 
contains 12 sites and the direct minimization of (1) for 
small samples of size multiple of 12 indeed reveals a Neel 
long-range order with 12 noncoplanar sublattices point- 
ing toward the 12 centers of edges of a cube (Fig. 2). The 
apparent complexity of this structure is somewhat light- 
ened when one considers the six spins around an hexagon: 
they lie in the same plane and make an angle 7r/3 with 



B. Discrete symmetry breaking 

It is clear that 0(3) is fully broken in the cuboc phase 
at T = 0: the point group symmetry of a cuboctahedron 
is simply that of the cube, i.e., Oh = Ox{Id,i} , where 
O is the octahedral group containing the 24 rotations 
leaving a cube or an octahedron invariant, and Id and i 
are, respectively, the identity and the spin inversion, or 
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spin flip. 

An important result arises when one considers the ac- 
tion of the spin flip alone. It is clear that i acting on a 
cuboctahedron takes it onto another cuboctahedron but 
the labels of the two cuboctahedra cannot be made to 
coincide by means of a global rotation. Namely, the or- 
der parameter we obtain is the mirror-symmetry image 
of the previous one (Fig. 4) . The order parameter in the 






FIG. 4: (Color online) An order parameter and its image by 
a spin flip. Clearly, after a global rotation of the sublattice 
directions, the two order parameters are related one to the 
other by a mirror symmetry. 

cuboc phase thus breaks the spin flip symmetry, and we 
are able to divide the ground states manifold into two 
classes. This makes us expect a transition at finite tem- 
perature associated with this Z 2 symmetry breaking. 
To show that this is indeed the case, we define a vari- 
able labeling the two classes of ground states. Consider 
the normalized scalar chirality on a triangle, namely 
■gj-r with (i,j,k) labeling the three sites 



(SiAS 

a A = 



clockwise. An inspection of the order parameter of Fig. 2 
reveals that oa is alternatively +1 on upward triangles 
and -1 on downward triangles (Fig. 5). We naturally de- 
fine the alternate scalar chirality as 



A 



1) 



^A, 



(2) 



where the sum runs over the 2N/3 triangles of the kagome 
lattice, and oa is, respectively, and 1 on upward and 
downward triangles. The spin flip trivially changes the 
sign of cr a so that m a — ±1 at zero temperature, de- 
pending on the class of the order parameter. Hence, mo- 
is the order parameter associated with the spin flip sym- 
metry breaking. 

Monte Carlo simulations have been performed on sam- 
ples of up to 1200 spins: they show that m CT vanishes at 
finite temperature while the associated chiral susceptibil- 
ity, defined as 



2N < 2 
k BXa = 7^r(< m a > 



< \m„\ > 2 ), 



(3) 



and the specific heat, seemingly both diverge (Fig. 6). 
Typical simulations involved 10 6 Monte Carlo steps per 
spin. The first results indicate that the transition is not 




FIG. 5: The normalized scalar chirality cta computed on the 
order parameter of Fig. 2. The spin flip trivially permutes the 
+ and — . 



in the two-dimensional Ising universality class. The com- 
plexity of the global set of excitations that does not re- 
duce to those of the above-mentionned oa variable are 
probably at the origin of a more complex behavior, pre- 
sumably a weak first-order phase transition 18 . The com- 
plete study of this classical phase transition will be pub- 
lished elsewhere 19 . 
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FIG. 6: (Color online) Specific heat C/ks (red dashed line) 
and chiral susceptibility ksXcr (blue solid line) versus temper- 
ature on a 1200 spins sample for J2/I Ji | — 0.5: both quantities 
diverge at T c = 0.24, simultaneously with the vanishing of the 
chiral order parameter m a (not shown). 



III. SPIN-1/2 MODEL 

We now turn to the spin-1/2 quantum model. Con- 
sidering the classical analysis, the question is whether 
quantum fluctuations are strong enough to wipe out the 
classical order. If so, we end up with a purely quantum 
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phase such as a valence bond solid or a RVB liquid 5 . If 
not, this means that quantum fluctuations merely dress 
the classical order parameter, reducing the mean sublat- 
tice magnetization while preserving its symmetries. Such 
states may thus be referred to as semiclassical states. 



A. A la Neel SU(2) symmetry breaking 



of the two groups. Indeed, if we restrict ourselves to the 
rotational symmetry breaking, we see that the original 
50(3) symmetry of (1) is reduced to its subgroup O in 
the Neel-ordered ground state. 

Thus, while D$ is an irreducible representation (IR) of 
50(3) of spin 5, it is an a priori reducible representation 
of O that one can decompose onto the five IRs T v of O 
according to 



The exact spectra of finite size samples have already 
been shown to be a powerful tool to find eventual Neel or- 
ders in quantum spin systems 20 ' 21 . Of course they allow 
for the exact computation of the square of the sublat- 
tice magnetization or any relevant structure factor. But 
mostly, while the extensive use of the symmetries of the 
Hamiltonian is compulsory in order to diagonalize large 
enough samples (typically 36 spins today) , it also gives a 
clear signature of SU (2)-breaking phases, even on small 
samples spectra. 

Namely, for an SU(2) symmetry breaking we expect a 
large set of low- lying eigenstates of the Hamiltonian, with 
different total spin 5 values, to collapse onto the ground 
state when the size of the sample N — > oo. These states 
have been called the quasidegenerate joint states (QDJS) 
and they have been recently observed experimentally 
on nanomagnets 22 . If the semiclassical picture is valid, 
these QDJS are expected to have energies well below the 
magnon excitations and to scale as 5(5 + l)/iV, as ex- 
pected for a quantum top. Exact spectra are thus dis- 
played versus 5(5 + 1) and the QDJS are often referred 
to as the Anderson's tower of states 21 . To have a true 
SU(2) symmetry breaking in the thermodynamic limit, 
the QDJS should have total spins up to 5 ~ 0(VN). The 
breaking of the SU(2) symmetry then occurs with all 
the QDJS collapsing onto the absolute ground state like 
1/N when N — > oo, i.e., faster than the softest magnons, 
whose energies scale as l/\/N, defining the a la Neel 
5£/(2)-symmetry breaking scheme 21 . This result explains 
why it is numerically more favorable to look at the QDJS, 
since order parameters only scale as 1/ y/N. 
In the thermodynamic limit, the ground state is a super- 
position of an infinite number of QDJS with different 5 
values, which clearly breaks SU (2) . 
The crucial point is that in each spin sector the number 
and symmetries of the QDJS are exacly determined from 
group representation theory 3 ' 21 by the symmetry of the 
expected ground state. 



B. Determination of the QDJS 

If the thermodynamic ground state exhibits a semi- 
classical 12-sublattice Neel order then the QDJS, if ever 
they exist, should be of symmetry compatible with both 
those of the Hamiltonian (1), since they are eigenstates, 
and those of the cuboc phase, i.e., Oh ■ 
A classical result of group theory is that the number of 
such states is completely determined by the structures 



v=l 



with 



(4) 



(5) 



geo 



where x v (g) an d Xs(flO are the characters of g € O in 
the IR T v of O , and in the IR D s of SO (3) , respec- 
tively. The character table of O is given in Table I and 
Xsid) ~ Sln ^n(e/l)^ 2 ^ ' w ith 9 the angle of the rotation g. 
For completeness we give the explicit decomposition (4) 
for spins up to 5=6 in Table II. 



TABLE I: Character table of O . Its irreducible represen- 
tations are named both with the usual nomenclature and 
with the numbering 1 < v < 5 used throughout the article. 
Also shown is the character table of the group {Id, i} . The 
character table of the whole Oh group is obtained by forming 
the direct product of the two tables. 
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Id 8C 3 3C* 2 6C7 2 6C* 4 
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1 1 


1 


-1 -1 
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2 -1 


2 
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Ti 
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-1 


-1 1 
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T 2 


3 


-1 


1 -1 


5 



{Id,i} 


id i 


r e 
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i i 
i -i 



The decomposition (4) directly gives the number of 
states that belong both to Ds and T v , i.e., that are 
compatible with both the SO{3) and O symmetries, 
as required for the QDJS. It should be emphasized 
that for a given S value the number of such states 

5 



is n(S) 



n v {S) dimr„ =25+1, as expected for a 



complete 50(3) breaking. 21 

We stress that (4) relies only on group theory and that 
it makes no reference to the representation space 24 . 
Now, if we are to find the total content of each spin sec- 
tor of the Anderson's tower of states, we should treat the 
whole Oh group, not limiting ourselves to the rotational 
symmetry-breaking 50(3) — > O as in (4). 
This is particularly simple since Oh is the direct product 
of O with the group {Id, i} , and the spin flip being also 
a symmetry of the Hamiltonian, there is no compatibility 
issue here. Thus, everytime an IR T v appears in (4) we 
actually get two copies of it associated to the two IRs 
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TABLE II: The decomposition (4) explicited for S < 6. 



D 


= A l 
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= Ti 
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T 2 




D 4 


= A 1 
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Ti 


+ T 2 


D 5 


= E 
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2Ti 
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T 2 




D 6 


= A 1 


+ 


A 2 


+ 


E 


+ Ti + 2T 2 



of the {Id, i} group (Table I). Since these two IRs differ 
only in their parity under the spin flip operation, which 
itself transforms any order parameter into its Z 2 image, 
this double quasidegeneracy is clearly reminiscent of the 
Z 2 symmetry breaking observed classically. 
We thus have formally determined the number and sym- 
metries of the QD JS appearing in each spin sector of the 
tower of states. However, for a given total spin, their sym- 
metries are given in terms of the IRs of Oh while exact 
diagonalizations provide eigenstates of given symmetry 
under the lattice symmetry group. 

It thus remains to map the IRs of Oh onto those of the 
lattice symmetry group, namely Gn — 2~W A Pn, where 
Tjy contains the N/3 translations by a Bravais lattice 
vector and Pn is the point group of the sample (in gen- 
eral, Pjv is a subgroup of Cq v , the point group of the 
infinite kagome lattice). 

Such a mapping clearly exists since the labeling of the 12 
vertices of the cuboctahedron induces a labeling of the 
lattice (Fig. 2). Applying an element of Oh to a cuboc- 
tahedron means permuting its 12 labels, which itself is 
equivalent to a lattice transformation. 
We thus have a mapping between Oh and some elements 
of Gat . Note that while the mapping between the group 
elements is not necessarily one-to-one, and in fact it is 
a one-to-many mapping for all the elements of the sub- 
group O , the resulting mapping between the IRs of the 
two groups is one-to-one, as is explicited in Table III. 
A notable exception is the spin flip, which is exactly 
mapped onto the rotation of the lattice R„ by angle tt 
around the center of an empty hexagon. Hence, the par- 
ity of an eigenstatc of the Hamiltonian (1) under is 
directly equal to its parity under the spin flip. Thus, if 
ever the thermodynamic ground state has the classical 
12-sublattice structure, we are able to find the number 
and degeneracies of the expected QDJS using Table II 
and the one-to-one mapping between the IRs of Oh and 
those of G N (Table III). 

However, there is still one subtlety. Indeed, in order not 
to artificially frustrate the 12-sublattice order, we choose 
samples containing multiples of 12 spins, i.e., N=12, 24, 
and 36 spins. The representation spaces of these three 
samples have different properties since the total spin on 
each sublattice is N/24, which may be integer or half- 
integer. 

To be more specific, we want to write down the matrix 
U (g) G SU(2) associated with a particular rotation g&0 



TABLE III: One-to-one mapping between the IRs of Oh 
and those of Gn ■ We considered a sample having the three 
edge centers in its first Brillouin zone and the full Ce v point 
symmetry. Ca v is generated by -R 27r /3, the rotation of the 
lattice by angle 2tt/3 around an empty hexagon, and a, a 
reflection whose axis is a diameter of an empty hexagon, so 
that the IRs of Ce v may be labeled by the associated three 
quantum numbers TZ 2n / 3 , Kjtj and a. The IRs of TV are 
labeled by the k vector of the first Brillouin zone. 



Q \ k TZ 2 n/3 gjV 
Ai 1 11 
A 2 1-12 

E j,j 2 3 
Ti Xi,2,3 -1 4 

T 2 Xi,2,3 1 5 



that acts on the wave function \cuboc) of a cit&oc-ordered 
state. The Hilbert space that contains such states is a 

12 

subspace of at/24, where T> N /24 is the Hilbert space 
i=i 

of one spin N/24. Thus, a natural choice for U(g) would 
be the tensor product of 12 Un /24(g) matrices, each one 
of them representing g in T> N / 24 . 

Now, we know that in each subspace D N / 2 4, if N/24 is a 

half-integer, Un /24(g) an d ^U^/ 2 4(g) are equally suitable 
choices, and we cannot decide between them other than 
arbitrarily, due to the double connectedness of 50(3) . 
Let us choose 12 such matrices anyway and form their 
tensor product U(g). When acting on a particular order 
parameter, named \cuboc) , we get 

U(g)\cuboc) = ip(g)\cubod), (6) 

where \cuboc') represents the order parameter obtained 
by applying the global rotation g^O on \cuboc) , and 
<p(g) is an overall phase factor that we cannot get rid of 
since it embeds the arbitrariness of our choice of the ma- 
trices U N /24(g) ■ 

Recall, however, that U(g) acts in a subspace of 
12 

T> N /24, which is known from spin algebra to contain 

i=i 

states with integer spins only, whatever N/24, integer, or 
half-integer. Hence U(g) is always a true representation 
of g, i.e., no ambiguity should remain in it and conse- 
quently in <p(g). Thus the group law should be exactly 
verified by ip(g), which is then simply a one-dimensional, 
thus irreducible, representation of O . 
Direct computation for N=12, 24, and 36 spins indeed 
shows that ip(g) — Xv (g)i where vq=1 for A^=24 spins, 
and i/o=2 for N=12 and 36 spins. 

Finally, to embed this subtlety in the decomposition (4), 
one simply has to permute A\ <-> A 2 and T\ <-» T 2 in Ta- 
ble II for N—12 and 36 spins, as can be seen directly 
on the character table of O (Table I). Note that no ad- 
ditionnal ambiguity arises from the consideration of the 
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{Id, i} part of Oh , since the whole discussion relies on 
the double connectedness of 50(3) that is irrelevant here. 

C. Analysis of exact N=36 spectrum 

We determine the number and degeneracies of the 
QDJS expected for a sample with 7V=36 spins and com- 
pare this to the exact diagonalization result. 
We first compute the decomposition (4) for total spin 
5 < 6, the approximate maximum spin of the QDJS for 
N = 36. As stated earlier, we just have to take Table II 
and perform the relevant permutations A\ <-> A 2 and 
T\ <-> T 2 . In order to take the whole Oh group into ac- 
count we recall that each one of the IRs of the decompo- 
sition should actually appear twice with the two possible 
parities under the spin flip operation. Then we map these 
IRs onto those of Gn and here we need to specify the 
shape of the sample we diagonalize (Fig. 7). The map- 




FIG. 7: A=36 sample. 



ping is readily done since the iV=36 sample has the full 
Cq v symmetry and we may directly use Table III. Hence 
we have obtained the full composition of the Anderson's 
tower of states in the lowest spin sectors < 5 < 6 (Ta- 
ble IV). Again we stress that each QDJS appearing in 
Table IV should be present twice, with the two parities 
under the operation, this last result being the hall- 
mark of the Z 2 symmetry breaking observed classically. 



The comparison with the exact spectrum of the Hamilto- 
nian (1) for the ./V=36 sample is straightforward (Fig. 8 
and TABLE V). Indeed, one clearly notes a set of low- 
energy eigenstates, well separated from the magnon ex- 
citations for total spin 5 <6, which scale reasonably as 
5(5 + 1). But the main point is that for each total spin 5 
the number and symmetries of the low-lying eigenstates 
are exactly those obtained from our symmetry analysis, 
as can be readily verified in TABLE V. 



In particular, we note that we have 2(25 + 1) QDJS in 
each total spin 5 sector, consistently with a complete 



TABLE IV: Expected symmetries (fc: wave vector, Tl^/i- 
phase factor in a lattice 2-k/3 rotation, and a: phase factor 
in a lattice reflection) and dimensionality (d) of the IRs 
appearing in the QDJS of the 12-sublattice Neel state. Each 
of these IRs actually appears twice in the spectra with 
TZtt — ±1 (not shown), these two copies will be noted in the 
following e (for even) and o (for odd). The last columns give 
the number of each IR expected in the S sector (up to 5=6) 
for the A=36 sample, according to Eq. (5). The last line 
gives the number n(5) of QDJS in each sector [counting the 
Z2 replica gives a total number of 2 n(S) QDJS in the spin S 
sector] . 
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1 3 5 7 9 11 13 



TABLE V: Energies per spin and symmetries of the lowest 
eigenstates of the Hamiltonian (1) for 5 < 4. 



s 


E/N\,h\ 


IR 





-0.420920 


2e 





-0.420054 


2o 


1 


-0.418935 


5o 


1 


-0.418596 


5e 


2 


-0.415654 


3e 


2 


-0.415317 


4o 


2 


-0.415282 


4e 


2 


-0.414885 


3o 


3 


-0.410506 


le 


3 


-0.410423 


4o 


3 


-0.410403 


5o 


3 


-0.409898 


4e 


3 


-0.409886 


5e 


3 


-0.409738 


lo 


4 


-0.404181 


3e 


4 


-0.403855 


5e 


4 


-0.403797 


2e 


4 


-0.403574 


4o 


4 


-0.403272 


4e 


4 


-0.403172 


5o 


4 


-0.402928 


2o 


4 


-0.402888 


3o 



SU(2) breaking in the thermodynamic limit, with the 
factor 2 coming from the two replica TZ n = ±1 of each 
QDJS and taking care of the Z 2 symmetry breaking. 
The same analysis has been made for N=12 and 24 non- 
frustrating samples, leading to the same result. 
These symmetry arguments are strong enough to claim 
that, at least in a certain range of parameters, quantum 
fluctuations do not destroy the complicated 12-sublatticc 
classical long-range order and that there exists a quan- 
tum cuboc phase in the thermodynamic limit, in the sense 
of a ground state consisting of the classical cuboc state 
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FIG. 8: Top: Exact spectrum of the Hamiltonian (1), with 
J2/IJ1 1=0.5, for the JV=36 sample shown in Fig. 7. The exact 
energies per spin are displayed versus S(S + 1). The + (resp. 
— ) symbols are eigenstates with even (resp. odd) parity un- 
der the TZn lattice rotation. Other symbols are eigenstates 
without the IZ^ symmetry. The full symmetries of the low- 
est eigenstates are given in Table V. The expected QDJS for 
the 12-sublattice Neel order appear between the two dashed 
lines. In particular, even and odd parity eigenstates come in 
an equal number in each spin sector, as expected for the Z2 
symmetry breaking. Bottom: Zoom on the QDJS (energies 
rescaled) . 



state. Thus, in this frame the classical cuboc ground state 
is simply a ferromagnetic state to which one can readily 
apply the Holstcin-Primakov transformation. 
First, we choose a particular cuboc state, say, the one of 
Fig. 2. At a given site i of the kagome lattice we define 
the local frame (xi, yi, Zi) with Zj the unit vector par- 
allel to the local classical spin S,. To choose Xi we note 
that each site of the kagome lattice belongs to two tri- 
angles pointing toward opposite directions. Consider the 
other two spins on the downward triangle and label them 
j and k, with (i,j,k) turning clockwise. The directions of 
the two spins j,k in our cuboc state are zj and z^ and 
one can easily verify that Xi = (z/j — Zj)/V2 is indeed a 
unit vector orthogonal to Z,. We completely determine 
the local frame by imposing y^ = Zi A Xi. 
This construction is translationally invariant and repeat- 
ing it for all the sites of the kagome lattice will lead 
to 12 different local frames associated to the 12 sublat- 
tices of the classical ground state. Hence, using the ap- 
propriate transition matrix IZi from the reference frame 
(x,y,z) to the local frame (xi,yi,Zi), we may compute 
the components of the spin at site i in its local frame 
St = (Sf ,Sf ,Sf ) from S { = TZ^ (Remark: there are 
only 12 different IZi matrices). 

Before computing the Hamiltonian (1) in the local frame, 
we note that it can be rewritten as a sum over the iV/3 
empty hexagons (Fig. 9): 



h = Ji y^ y^ Si ■ Sj + J2 y^ y^ s^-s* 



(7) 



O <i,j> 



O <<i,k» 



where <i,j> and «i,k» are now, respectively, 
the six nearest and six next-nearest neighbor pairs of 
sites enclosed in the empty hexagon O. Now, using 
Si ■ Sj — S'jTijS'j, with Tij = t TZ{R,j 1 we may rewrite (7) 
in the local frame as 



renormalized by quantum fluctuations, as explained ear- 
lier. 

However, it can be objected that long-wavelength quan- 
tum fluctuations, which cannot be accounted for on the 
small samples we diagonalized, may wipe out the cuboc 
order. 



IV. SEMICLASSICAL APPROACH 



O <i,j> 



SiTijSj 



E 

O <<i,fe>> 



S'iTikS'k- (8) 



Now, we are able to quantize the fluctuations around the 
classical cuboc state using Holstcin-Primakov bosons. Up 
to quadratic order the transformation at site i is written 
as 



S'+ = S^+jSf* = V2S--n,r; 
S'r = S^-jSr = 4V2S~- 



~ V2Scl (9) 



m 



Now that we are convinced that the 12-sublattice Neel 
order observed in the classical cuboc phase also exists in 
small samples of spins 1/2, we may compute the effect 
of long-wavelength quantum fluctuations on the energy, 
sublattice magnetization, and chiral order parameter in 
the ground state using the spin-wave approximation. 
The route to compute the quantum deviations to a clas- 
sical cuboc state is straightforward: at each site of the 
kagome lattice wc define a local frame in spin space whose 
z axis is aligned with the local spin in the classical ground 



where c\ and Ci, respectively, create and annihilate a 
Holstein-Primakov boson at site i, with S the length of 
the classical local spin and tii — cjcj. 
Inserting (9) into (8) we obtain the quantized version of 
the original Hamiltonian (1) up to quadratic order. As 
usual we Fourier-transform a and c\ according to 



1^' 



(10) 
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FIG. 9: (Color online) Kagome lattice with nearest (solid) 
and next-nearest neighbor exchange (dashed). Note that we 
span all the exchange paths by considering the 12 links inside 
an empty hexagon o located at Rq and then iterate over the 
whole Bravais lattice, u and v are the basis vectors of the Bra- 
vais lattice of length 2 and the three sites a, f3, 7 per Bravais 
cell define the three unit vectors e a — u/2, ep — (v — u)/2 



and e. 



-v/2. 



where the sum runs over the first Brillouin zone, Rq is 
a vector of the Bravais lattice, and /ij = a, 0, 7 indicates 
one of the three possible sites in the Bravais cell (Fig. 9). 
Again, in the local frame the complicated 12 sublattice 
order is just a ferromagnetic state, so that we need only 
three flavors of bosons fii = a, /3, 7 associated with the 
three sites per Bravais cell on the kagome lattice. We 
may finally bring the Hamiltonian to matrix form as 



n = (Ji- j 2 )ns{s + 1) + E V}M q V q , 



(11) 



where the sum runs over the entire first Brillouin zone, V q 
is the column vector (c^, c^, c^, c_g, c?} q , c_„), and M q 

is the 6 x 6 matrix M q = (J 2 - J X )S1 + ( „ 9 f q 

\ a -q A -q 

with 1 the identity matrix, and 

a a p{q) a aJ (q) 
A q = \ a aP (q) a Pl {q) | , (12) 
a aj (q) a^{q) 



and 



b a p{q) b aj (q) 
B q = \ b a0 (-q) b Pl (q) } , (13) 
b ai (-q) b 07 (-q) 



with 



, . J2S JiS J\S . 

a a /3{q) = —j— cos <2-/3_ Q — cosg 7 - —j= smg 7 , 

. . J2S JiS J2S . 

b a p{q) = —^-cosqp-a — cosg 7 + — —smq - a , 

with the four remaining matrix elements obtained by 
cyclic permutation of the indices and where we have used 
the condensed notation qp- a = q ■ (ep — e a ). 



Note that the first term in (11) contains the usual dom- 
inant contribution to the renormalization of the classi- 
cal ground state energy (Ji — J 2 )NS 2 . Note also that 
A q = A q and B q — B_ q so that M q is indeed hermitic. 
As usual, the next step is to find a transition matrix P q 
such as M q is diagonal in the new basis. 
It should be emphasized that P q is strongly prescribed 
by the fact that it has to preserve the boson commuta- 
tion relations 23 , much in the same way as the Bogoliubov 
tranformation for the collinear antiferromagnet on the 
square lattice. 

We end up with three Bogoliubov bosons, obtained 
from Wq={d a q , dP r d? q , d% d^, (fJ q )=P q V q , and the 
six eigenvalues give the corresponding three dispersion 
branches. 

In the new basis (11) reads as 

n = (j 1 -j 2 )NS(s + i) + J2 E <«< + 5)> 

q fi=a,f3,y 

(14) 

so that the energy per spin in the ground state, which is 
the vacuum |0) of Bogoliubov bosons, is simply 

e» = 1{Q\H\0) = (J 1 -J 2 )S(S+1) + ±Y1 E % 



(15) 

As for the collinear antiferromagnet, note that the Bo- 
goliubov transformation P q is actually singular wherever 
lu^ = 0, namely at q ~ X\^ 2 ^, where soft modes are ex- 
pected in the thermodynamic limit. One indeed verifies 
that the lowest branch vanishes at each one of the three 
edge centers of the first Brillouin zone giving exactly 
three Goldstone modes in the thermodynamic limit, as 
expected for a complete SU(2) symmetry breaking. 
We may then compute the renormalization of the mag- 
netization in the local basis, 



1 



JV 



= ^5< lE 5 f|0> 

i— 1 



^a-^E'E^n 



. (16) 

where the prime denotes a sum over the first Brillouin 
zone deprived of its three edge centers and where P^ is 
the (i, j)th matrix element of the matrix P q . 
Another quantity of interest to us is the renormalization 
of the scalar chirality on a triangle. It is naturally nor- 
malized by its value in the classical ground state, so that 
we define £a = -0(S*i A Sj) ■ Sk on each triangle (i,j,k), 
with (i,j,k) turning clockwise. We may then compute the 
renormalization of the alternate scalar chirality in the 
ground state as 



,N 



-4<°iE(- 1 ) aA ^io) = i + |(i-^E'^)' 

A q 

(17) 

where the sum runs over the 2N/3 triangles of the kagome 
lattice, while cva is, respectively, and 1 on upward and 



8 



downward triangles, and 

3 



( cos q a (P^ + P^)(P^ +P$' j )\ 



cos^(P^' + P^')(P^ + P, 6 



cos q 7 (P^+P^)(P t 



2,3 



p5,j) 
1 > 



-y- sin g c 



(p5,jp6,j _ p2,jp3,j\ 
\ q q q q J 



^ sin q0 (p*JpSj -pldpK 



V+ # sin g 7 (P^ P q 5 



5,j 



p^ p^ ) y 



These quantities were numerically computed on finite 
size samples with linear sizes L < 10 2 Bravais lattice 
spacings. As usual, the leading correction to the classi- 
cal value comes from the first magnon excitation whose 
energy scales as 1/L. We find very good agreement with 
the expected scaling laws for (~ l/£ 3 ), m N (~ 1/L), 
and (~ l/£), and perform the extrapolation to the 
thermodynamic limit. 

We first note in (Fig. 10) that in the thermodynamic limit 
m°° and m|° remain finite in a large region around the 
point were exact diagonalizations were performed, thus 
showing the stability of the cuboc phase against long- 
wavelength quantum fluctuations. For J2/I Ji| = 0.5 we 
find to lowest order in the spin wave approximation that 
m°° and are renormalized by 16% and 50%, respec- 
tively. 

However, upon increasing J2, both quantities decrease 
drastically. For J2/IJ1I > 3 the chirality disappears 
which casts a strong doubt on the stability of the 12- 
sublattice state itself in this range of parameters. 




FIG. 10: (Color online) Extrapolated values of the magne- 
tization in the local basis m°° (solid line) and the alternate 
scalar chirality (dashed line) as a function of J2/IJ1I. 



-2.05 




V. A GAPPED PHASE FOR J 2 /\Ji \ = 5.0 

Exact spectra for J2/IJ1I = 5.0 differ notably from 
spectra of an ordered phase as described in Sec. III. First 
of all, their low-lying levels in each 5* sector do not scale 
as S(S + 1) but rather as S, as can be seen in Fig. 11. 
Second, the gap to the first excitation seemingly does not 
close to zero with the system size, as shown in Fig. 12. 
This second result is consistent with the first one: in 
a gapped phase a finite magnetic field H c is needed to 
close the gap. This critical field is directly proportional to 
the first derivative of the energy versus S. For this given 
model in this range of parameters, the available sizes of 
exact spectra are large enough to infirm the presence of 
a 12 sublattice Neel order and confirm the existence of a 
gapped phase. However they are too small to discriminate 
between a true spin liquid or a valence bond crystal, and 
to decide if there is one or two different gapped phases 
in this range of parameter. 



FIG. 11: An exact spectrum of the iV=36 sample for 
^2/1^/1 1=5.0 versus total spin S. The symbol convention is 
the same as that of Fig. 8. 



VI. CONCLUSION 

In this paper we have studied the Ji — J 2 model on the 
kagome lattice with J\ < and J 2 > 0. We have found a 
12-sublattice ordered phase for J\ < and J2/IJ1I > 1/3. 
This new phase was shown to resist quantum fluctua- 
tions. On the exact spectra of small size samples we found 
the complete signature of this complicated Neel order, 
i.e., the number and symmetries of the QDJS in the tower 
of states, based on a very general group-theoretical ap- 
proach. Moreover, in the spin wave approximation, long- 
wavelength quantum fluctuations were found to renor- 
malize the order parameter to a finite value in a finite 
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FIG. 12: (Color online) Left: Finite size scaling of the spin 
gap AE for «7a / 1 «7i | = 5.0, in the supposed to be gapped 
phase. Right: The same gap in the 12-sublattice Neel phase, 
for J2/I Ji\ = 0.5, closes as 1/N, as expected for the a la Neel 
SU(2) symmetry breaking (Sec. Ill A). 



range of parameters up to J2/I J\ \ — 3. 



The noncoplanarity of the 12 sublattices in the cuboc 
phase was shown to induce a chiral symmetry breaking, 
to which we associated a chiral order parameter. Classi- 
cally, we showed that this Z2 symmetry was restored at 
finite temperature, consistently with the Mcrmin- Wagner 
theorem, though the exact nature of the transition re- 
mains to be investigated. We were also able to find the 
signature of this discrete symmetry breaking on exact 
spectra. 

Finally, we showed in the spin- wave approach that the 12 
sublattice order is wiped out by quantum fluctuations for 
^2/l^rl > 3. Exact diagonalizations indeed confirm the 
existence of a spin-gap phase, with short-range order in 
spin-spin correlations. For this model the largest available 
sizes (TV = 36) are too small to give more information on 
the nature of this quantum phase. 

The Laboratoire de Physique Theorique de la Matiere 
Condensee is UMR 7600 of the CNRS. 
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